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Abstract 

The paper is devoted to the calculation of the two-loop 0(a B y 2 ) and 0(y 4 ) MSSM 
corrections to the relation between the pole mass of the 6-quark and its running mass in 
DR scheme. To evaluate the needed diagrams the large mass expansion procedure is used. 
The obtained contributions are negative in most of the regions of the parameter space and 
partly compensate the positive (D(a 2 ) contribution calculated earlier. 

1 Introduction 

The Minimal Supersymmetric Standard Model (MSSM) became an underlying framework for 
both the theoretical and experimental research in supersymmetric phenomenology. During re- 
cent years considerable progress has been achieved in calculation of various processes involving 
supersymmetric particles and the parameter space of the MSSM is severely constrained by nu- 
merous phenomenological fits ^ El EH El E] • These fits require precision calculations of radiative 
corrections due to virtual supersymmetric particles. Because of complexity of the MSSM La- 
grangian these calculations are very lengthy and assume the creation of the proper computer 
codes. 

In this paper we calculate two-loop radiative corrections to the pole mass of the 6-quark 
within the MSSM. The one-loop correction has been calculated in |B] and is very essential. It 
reduces considerably the value of the running quark mass extracted from the pole mass. It 
should be noted that the leading 0(a s ) contribution from stop-gluino loop which is positive is 
reduced significantly by large negative correction from stop-chargino loop that is proportional 
to the Yukawa couplings of heavy quarks (ytiUb), tan/3 and the Higgs mixing parameter [i. For 
large tan/3 one has to take into account not only y t but also yb since it becomes large as well. 
Calculation of the two-loop correction of the order of 0(a 2 s ) has been performed in 0. It was 
found that this correction is also positive and of the same order of magnitude as the one-loop 
MSSM contribution. 

The aim of this paper is to present our calculations of the corrections proportional to the 
Yukawa couplings. We calculate the two-loop 0(a s y 2 ) and 0(y 4 ) MSSM corrections to the 
relation between the pole mass of the 6-quark and its running mass in DR scheme jS] . The DR 
scheme corresponds to a theory regularized by dimensional reduction (from space-time dimension 
d = 4torf = 4 — 2s) and renormalizcd minimally. This allows one to keep supersymmetry 
unbroken and use the computational advantages of dimensional regularization. 

To evaluate the needed diagrams we made use of the large mass expansion procedure |S] 
and restrict ourselves to the terms up to 0(m 2 / M 2 ard ), where Mhard stands for all mass scales 
involved in the problem that are much larger than mb- Since the total number of diagrams 
exceeds 1000, we are not able to present the calculation in a compact form, so we provide 
numerical demonstration of the results. The final result is also available upon request in a form 
of C++ code 1 . As expected the calculated corrections are negative in most of the regions of the 

1 E-mail: mailto: vargOthsunl .jinr.ru. 



1 



parameter space of the MSSM and partly compensate those proportional to 0(a 2 s ). By absolute 
value the "leading" 0(a 2 s + a s y 2 + y 4 ) two-loop corrections are of the order of 30 to 40 percent 
of the one-loop ones. 



2 Pole mass of 6-quark 

The pole mass of a particle is defined as the real part of complex pole of resumed propagator 
(we discuss only perturbative effects). Full connected propagator of a quark can be written as 2 

(1) 



p — TO — TOi) ' 

where 

S(p, mi) = pT< v (p 2 ,m 2 ) +p7 5 XU(p 2 ,TO 2 ) + m^sip 2 ,mi) (2) 
is the self-energy of the quark, so the pole mass M po i e satisfies the following equation 

(1 + Z v (M 2 ole , to 2 )) 2 M 2 olE - Y? A (M 2 pole , m 2 )M 2 ole - to 2 (l - E s (M 2 oie , to,)) 2 = 0. (3) 



Mpoi e - TO 



Solving this equation perturbatively, one gets 

= aM {V > +a 2 M {2 \ where (4) 

III 

M« = E^(to 2 ,to 2 ) + E^V 2 ,™*), (5) 

M ^ = 4 2 V,™ 2 ) + 4 2 V 2 ,^) + ^ V,™?) 

+M« ^W(to 2 ,to 2 ) + 2to 2 A (e£V >«*?) + sgV.mO) Ip 2 =- 2 ) > (6) 

and a stands for all couplings of the theory. 

Using eq. 10} , one calculate relation between pole and running masses of the 6-quark (de- 
pending on the prescription used for evaluation of E™, mass parameter to in Q may correspond 
to bare or renormalized mass, so one obtains relation between pole and bare masses or between 
pole and running masses). 

We use regularization by dimensional reduction DR, a modification of the conventional di- 
mensional regularization, originally proposed in [H] and the same renormalization prescription as 
in ■ Therefore we calculate 

Am b _ Ml ole - m r b un 



r ti a 



nib rrib 

where M^ ole stands for the pole mass of the 6-quark and to™" = m&(/2) corresponds to the 
running DR mass of the 6-quark at the scale ft. 

To evaluate this quantity, we have to calculate more than 1000 two-loop propagator type 
diagrams (see Fig.^). In order to simplify our calculations, we used the t'Hooft-Feynman gauge 
(£w — l>£z = 1)- We also neglected mixing in the chargino and neutralino sectors, so only 
higgsino states were taken into account. This implies that all higgsinos have the same mass 
which is equal to Higgs mixing parameter \x. 

2 The CKM matrix is supposed to be diagonal. 
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Figure 1: Feynman diagrams contributing to quark self-energies, Here q, q%, g,g denote quark, 
squark, gluon and gluino fields respectively, H. = {ho, Hq, Aq, Go, H + , G + } corresponds to various 
Higgs fields and \ — {xo,X + } are higgsino states. Two-loop diagrams of the order 0(a 2 s ) are 
not presented here (see, e.g jjj). 



Due to presence of many different mass scales, the problem of exact evaluation of two-loop 
diagrams is rather complex. Exploiting the fact that 

m b < m t ,m ho ,mH„,m H +,mq,m i o,rn ji +,rn Go ,rn G +, (8) 

where mj and nit are masses of bottom and top quarks, respectively, rrih , ttih , rn H + are masses 
of the Higgs bosons, m^o is the neutralino mass, is the chargino mass, rriq are masses of 
different squarks, mG 01 mg+ are masses of pseudo-goldstone bosons, we use the method of large 
mass expansion 9 to reduce the evaluation of multi-scale two-loop diagrams to the calculation 
of two-loop vacuum integrals and products of one-loop on-shell propagator-type diagrams and 
one-loop bubble integrals. We include only terms up to 0(ml/M^ ard ), where Mhard denotes any 
mass of the right hand side of ijSJ ■ 

A general Feynman diagram Ft which depends on the large masses Mi , Mi , . . . , small masses 
mi, m,2, ■ ■ ■ and small external momenta pi,p2, ■ ■ ■ can be expanded as follows 

F r (pi, . . . ,Mi, . . . ,mi, . . .) = ^F vh {p,m)M 1 {p 1 ,rn)F 1 {M,m,p~ l ), (9) 

7 

where operator AA 1 {p~ l ,m) performs Taylor expansion in small external (with respect to sub- 
graph 7) momenta and masses. The sum runs over all asymptotically irreducible subgraphs of 
original graph T. 
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r 7 = r 



Figure 2: The diagram has no cuts composed of two light lines, therefore the naive expansion 
can be used. 




Figure 3: The diagram has one cut composed of two light lines, so non-trivial subgraph 71 has 
to be taken into account. 



For calculation of a two-loop propagator diagram that does not have cuts composed of two 
light lines up to 0(p 2 /M 2 ) one can use naive expansion, when sum in runs over only the 
trivial subgraph 7 = T (see Fig. [2J. In case of a diagram that has such a cut a non-trivial 
subgraph 71 has to be taken into account (see Fig. 02). 

Two-loop vacuum integrals can be recursively reduced to a master-integral |1 1 j by integration 
by parts method [T2| . 

Our calculation can be performed in a semi-automatic way. First, FeynArts is used 
to generate the diagrams. Then the contribution of individual diagrams to the relation Q is 
evaluated by means of C++ program TwoLoop, based on GiNaC [T2|. TwoLoop JS] performs 
large mass expansion according to ©, uses recurrence relations of to reduce two-loop vacuum 
integrals to the master integral, and performs expansion in e. Results of this calculation were 
analytically cross-checked using FORM program based on ON-SHELL2 package |16|. 

After evaluation of two- loop diagrams one should substitute them into @ and perform renor- 
malization in order to express the correction J7J) in terms of running DR parameters. Counter- 
terms were calculated and checked in the same way as in jj], e.g. our renormalization constants 
for the gauge and Yukawa couplings reproduce well-known (3- functions [P7| . 

3 Numerical results 

The analytical result of our calculation is very complicated due to presence of a large number of 
masses, and no phenomenologically acceptable limit seems to exist. Therefore we present here 
the numerical analysis of our results. 

While a detailed scan over the more-than-hundred-dimensional parameter space of MSSM 
is clearly not practicable, even a sampling of the three- (four-) dimensional CMSSM parameter 
space of mo, rni/ 2 and Aq (tan/3) is beyond the present capabilities of phenomenological studies, 
especially when one tries to simulate experimental signatures of supersymmetric particles within 
a detector. For this reason, one often resorts to specific benchmark scenarios, i.e., one studies 
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only specific parameter points or at best samples the one-dimensional parameter space ^1 

Numerical values of running SUSY parameters at the Mz scale have been calculated as a 
function of CMSSM parameters with the program TwoLoop ^3 after interfacing it to a slightly 
modified version of SOFTSUSY 20 in the framework of mSUGRA supersymmetry breaking 
scenario. We only consider /i > and large values of tan/3, since small tan/3 and negative /it 
seem to be excluded by experimental data |2J. 

Fig. |H shows m 1/2 -dependence of 0{a 2 s ) [7] and 0(a s y 2 + y 4 ) MSSM contributions to J7J). 
For comparison, full one-loop MSSM correction was plotted. Numerical value of the 0(a 2 ) 
contribution is comparable with the full one-loop result, so neglecting other two-loop MSSM 
corrections leads to a contradiction with the perturbation theory. The 0(a s y 2 + y 4 ) contribution 
is negative in a wide range of the CMSSM parameter space and reduces 0(a 2 ) contribution by 
approximately 30 — 80%. As a consequence 0{a 2 + a s y 2 + y 4 ) two-loop correction appears to 
be less than 40% of the one-loop contribution and therefore can be used in phenomenological 
studies of CMSSM. 

Varying Aq in a wide range from —1000 GeV to 1000 GeV does not change the behaviour of 
0(a 2 + a s y 2 + J/ 4 ) correction significantly, as one can see from Fig. 0] To make this fact more 
clear, several curves corresponding to different values of Aq were plotted in Fig. [5j 

In Fig. we present dependence of different two-loop MSSM contributions on too. For 
mi/2 > 500 GeV the change of m Q from 400 GeV to 1000 GeV yields change of 0(a s y 2 + y 4 ) 
correction less than 10%. On the other hand, dependence of 0(a s y 2 + y 4 ) correction on too 
becomes essential for relatively small m 1 / 2 (m 1 / 2 ~ 200 GeV). 

It should be noted, that in considered range of CMSSM parameters space 0(a 2 + a 8 y 2 + y 4 ) 
correction is positive. 



4 Conclusion 

In this paper we presented the results of calculation of the two-loop corrections to the relation 
between pole and running massess of the &-quark, proportional to the Yukawa couplings of 
heavy quarks (yt,yb)- We provided a numerical analysis of the value of these corrections in 
different regions of the CMSSM parameter space. Corrections presented here can be used in a 
renormalization group analysis of the Yukawa coupling unification. It was also found that nib 
corrections significantly affect low energy phenomenology where the &-quark enters |21l . The 
analysis given in |22| showed that the neutralino relic density is very sensitive to the mass of the 
fe-quark for large tan/3, thus the result considered in this paper may have important implications 
to the dark matter searches. We are going to study these issues and calculate more 0(a 2 s ) terms 
in the large mass expansion of the relation between the top pole and DR masses, which requires 
some improvement of our code. 
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Figure 4: Different MSSM contributions to Am b (M z )/m b (M z ). Here m = 1000 GeV, left 
column corresponds to tan/3 = 35, right column corresponds to tan/3 = 50. 
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tan/3 = 35 tan/? = 50 
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Figure 5: A$ dependence of 0(a 2 s + y 4 + a s y 2 ) two-loop correction to Arrib(Mz)/mb(Mz). Here 
too = 1000 GeV. Varying A$ from —1000 GeV to 1000 GeV does not change behaviour of the 
correction significantly. 



m 1/2 = 200 GeV to 1/2 = 700 GeV 
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Figure 6: Different two-loop MSSM corrections to Am,b(M z) / rrib{M z) as a function of too- Here 
too = 1000 GeV, A = 0, tan/3 = 50. Dependence of 0(a s y 2 + y A ) correction on m is essential 
for mi/2 = 200 GeV. On the other hand, for TO1/2 = 700 GeV the change of too from 400 GeV 
to 1000 GeV yields change of 0(a s y 2 + y 4 ) correction less than 10%. 
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